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High-energy heavy ion collider experiments at RHIC and LHC have revealed that relativistic 
hydrodynamic models describe the hot and dense quark matter quantitatively. In this study, I 
develop a novel dissipative hydrodynamic model at finite baryon density to investigate the net 
baryon rapidity distribution. The results show that the distribution is widened in hydrodynamic 
evolution, which implies that the transparency of the collisions is effectively enhanced. This 
suggests that the kinetic energy loss for medium production at the initial stage could be larger. 
Furthermore, the net baryon distribution is found sensitive to baryon diffusion, implying that 
dissipative hydrodynamic modeling would be important for understanding the hot medium. 



Xth Quark Confinement and the Hadron Spectrum, 
October 8-12, 2012 

TUM Campus Garching, Munich, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Lieence. 



http://pos.sissa.it/ 



Dissipative hydrodynamic evolution of hot quark matter at finite baryon density 



Akihiko Monnai 



1. Introduction 

Quarks and gluons are considered to be deconfined from hadrons at high temperatures and to 
form the quark-gluon plasma (QGP) [p. It had been naively expected that the plasma is weakly 
coupled because of the asymptotic freedom of QCD. High-energy heavy ion collisions revealed, on 
the contrary, that the hot QCD matter would be strongly coupled near the quark-hadron crossover. 
The experimentally observed hadronic yields turned out to be described quantitatively by relativis- 
tic hydrodynamic models with small viscosity at the Relativistic Heavy Ion Collider (RHIC) and 
the Large Hadron Collider (LHC). Modern hydrodynamic analyses take account of viscosity and 
fluctuation to improve the model accuracy However, the net baryon number is neglected in 
most of those studies even though it would be important in understanding the experimental results 
such as particle-to-antiparticle ratios. This is possibly because viscous hydrodynamic equations of 
motion are complicated in the presence of finite charge density and also the boost invariant ansatz 
is often assumed, by which the net baryon number conserved at forward rapidity is omitted. 

Baryon stopping has been used to quantify the transparency of collisions [|3|, Q]. At early 
fixed-target experiments in the Alternating Gradient Synchrotron (AGS) and the Super Proton Syn- 
chrotron (SPS), the mean rapidity of the baryon number is largely reduced from that of the in- 
coming projectile, suggesting that the nuclei stay near the center of mass after the collision. The 
later collider experiments at RHIC and LHC, on the other hand, show that the collisions become 
increasingly transparent at high energies. This naively implies that only part of energy is used for 
the production of a hot and dense medium for the latter cases. In this study, I investigate the ef- 
fects of dissipative hydrodynamic medium interaction at the later stage on the net baryon rapidity 
distribution to find that more energy could be available for the hot medium j5|, q). 

Development of a finite-density dissipative hydrodynamic model would also be useful in the 
search of the critical point on the phase diagram of QCD because the first principle approaches 
have difficulty at non- vanishing chemical potential due to the sign problem. The applicability of the 
hydrodynamic picture at lower energies would need investigating through comparison of theoretical 
calculations and experimental data since the beam energy scan at RHIC show that large momentum 
azimuthal anisotropy would stay relatively large [^J. In such systems dissipative process, which 
accounts for the deviation from local equilibrium, would become important. The beam energy 
scans are also planned and proposed at Facility for Antiprotons and Ion Research, Nuclotron-based 
Ion Collider fAcility and Japan Proton Accelerator Research Complex. 

2. Finite-density dissipative hydrodynamics in relativistic systems 

Dissipative processes in relativistic systems are non-trivial since the relativistic extension of 
Navier-Stokes equation is subject to acausality and instablity. The equation is obtained by intro- 
ducing off-equilibrium correction to the entropy current up to the first order in terms of dissipative 
currents. On the other hand, one can consider the next order corrections that introduce relaxation 
effects to the theory. This second order theory is known to be causal and stable and is widely used 
for the analyses of dissipative phenomena in relativistic fluids. In this study I employ the Israel- 
Stewart formulation [||] extended for the systems with particle number changing processes as 
the foundation of the causal dissipative hydrodynamic model. 
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The flow tensor decomposition of the conserved quantities yields thermodynamic variables; 
the energy density eo = u^u v Tq V , the pressure Po = —A^ v Tq V /3 and the net baryon number rigo = 
UpNgQ. The dissipative currents are the bulk pressure n = — A )Iv dT llv /3, the shear stress tensor 
7Z^ V = [(A tl a A v ji +A^A v ,)/2-A^ v A a p/3}8T a P = 8T^ V \ the baryon dissipation current V = 
A^ v Vg . Here u^ 1 is the flow and A^ v = g^ v — u^u v is the projection operator. Energy-momentum 
and baryon number conservations read 



d^ v = u v [De + {e + P + 11)^^-71^^} 

+ (eo + PQ + U)Du^-V^(P Q + n)-K^ v Du v + A^ v VPK vp = 0, 
Dn J0 + n J0 V^ + V^Vj 1 - VfDu^ = 0. 



9Xb 



(2.1) 
(2.2) 



Here D = u^d^ and V 1 = A^ v d v reduce to the time and the space derivatives in the local rest frame. 
The entropy production can be written in the extended second order theory [Q] as 



-ux n -v^+7z„ v xr>o, 



(2.3) 



where the thermodynamic forces with the next-to-leading order corrections are 
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(2.6) 



Here T is the temperature and ju# is the baryon chemical potential. The semi-positive definiteness 
of the entropy production suggests IT = — ^nn^n, V* 1 = ft/Xy and n^ v = 2t]X 7[ where £nn, Ky and 
T] are the (bare) bulk viscosity, the baryon charge conductivity and the shear viscosity, respectively. 
The linear cross coefficients - the bulk-energy and the bulk-baryon cross viscosities £n§ e and £ns« s 
and the baryon-heat conductivity Kyw - can be either positive or negative. The linear response 
theory require that the coefficients satisfy Onsager reciprocal relations [1C, 11], a representation 
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of microscopic reversibility. Tn, Ty and x n are the relaxation times and #'s are the second order 
transport coefficients. 

Since the main focus is on time-evolution of the system in the rapidity direction as the net 
baryon number originates from the shattered nuclei which resides at forward rapidity, transverse 
dynamics is simply integrated out in the model. This would be justified for the net baryon distri- 
bution because the dependence of the net baryon yields on transverse geometry is experimentally 
shown to be very small @. 

The equation of state (EoS) and the transport coefficients are necessary to solve the hydrody- 
namic equations of motion. Latest results of lattice QCD calculations for the EoS in the vanishing 



chemical potential [13] and the quadratic baryon fluctuation [14] are used for constructing the 
finite-density EoS in Taylor expansion method. Due to difficulties in first principle calculations 
of the dynamical responses, the viscosity and the conductivities are chosen as v\ = s/4n [0]» 
£ = 5[(l/3) — cl\r\ []T6[ ||] and K v = cv(dji B /dn B ) T l /2n JT7p. Here s is the entropy density 
and cs is the sound velocity. They are based on Anti-de Sitter/conformal field theory correspon- 
dence except for the bulk viscosity, which is non-conformal and is estimated here with the method 
based on the non-equilibrium statistical operator method. Note that the scalar linear cross terms 
are merged into this effective bulk viscosity. The baryon-heat cross conductivity is parametrized 
as Kvw = c vw[ n BQ/ ' {^0 + Pq)](5tiT Ky) 1 / 2 from the argument based on matter-antimatter symmetry, 
semi -positive definiteness of the transport coefficient matrix and dimensionality cy and cyw are 
parameters. The relaxation times are set to x n = (2 — \n2)/2%T , Tn = 18 — (9 In 3 — V3k) /24kT 
and Ty = \rv2j2nT according to the string theoretical approach Jl7|]. The other second order trans- 
port coefficients are set to vanishing for the moment to see the qualitative responses of the system 
to each dissipative current more clearly. 

The initial conditions are estimated from the color glass theory, a description of gluon satu- 



ration. Energy density profile is employed from the gluon distribution in Nara adaptation [18, 19] 



of Monte-Carlo Kharzeev-Levin-Nardi model [2C, 21] by averaging it over the transverse area as 



eo(lb,T],s) = [ \ / S ama )dEj / TQdr] s where r\ s is the space-time rapidity. Here the gluon energy distri- 
bution is assumed not to be much modified before the initial time To, which is chosen as 1 fm/c. Net 
baryon density profile is obtained from the valence quark distribution in the color glass medium 



2, 23] as ubo(Tq, t] s ) = {I / S aiea )dN B _ B / Todrj s . The dissipative currents are set to vanishing at the 



initial time to better illustrate the effects of the off-equilibrium processes. 
3. Results 

Figure [j] shows the net baryon rapidity distributions of the color glass-based initial state and 
the ones after ideal, viscous and dissipative hydrodynamic evolution for Au-Au collisions in RHIC 
at y/SNN = 200 GeV and Pb-Pb collisions in LHC at ^/snn = 2.76 TeV. Here the viscous hydrody- 
namic results include shear and bulk viscous effects and the dissipative hydrodynamic ones include 
baryon diffusion effect in addition to the viscous ones. The parameter in the transport coefficients 
are set to cy = 1 and cyw = for the moment. Freeze-out is calculated in Cooper-Frye formula 



[24] with the freeze-out temperature Tf = 0.16 GeV. Note that the post-hydrodynamic evolution 
is not considered here because the modification would be small in dilute media and also the net 
baryon number does not change in hadronic decay. 
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Figure 1: Net baryon rapidity distributions of the color glass-based initial state (solid lines) and of ideal 
(dashed lines), viscous (dash-dotted lines) and dissipative (dotted lines) hydrodynamic results at (a) RHIC 
and (b) LHC. The data points are the scaled results of the BRAHMS Collaboration. 



One can see that the hydrodynamic evolution widens the net baryon distributions as the out- 
ward flow flux is induced by the pressure gradient and that the distribution after the hydrodynamic 
evolution exhibits overall agreement with the scaled results of the net proton distribution experi- 
mentally observed by the BRAHMS Collaboration [25] at RHIC. Effects of viscosity and diffusion 
are suggested to be more visible at RHIC than at LHC. Viscosity tends to lessen the hydrodynamic 
effect because the longitudinal pressure is reduced as dynamical responses to expansion. Diffusion 
also steepens the distribution since gradients in the chemical potential field are induced by the net 
baryon number at forward rapidity. It is note-worthy that this suggest the finite density dissipative 
processes, which is naively believed to be small, might have visible effects on observables. 

Baryon stopping is quantified with the mean rapidity loss (8y) = y p — (y) @], where y p is 
the rapidity of the incoming projectile and 

w= /V%lW dy /f^M iy . (3.1) 
Jo ay I Jo ay 

The rapidity losses of the initial conditions are (8y} = 2.61 at RHIC and 3.36 at LHC for the current 
parameter settings. It should be noted they are dependent on parameters in the color glass model 
and can be chosen differently. Ideal, viscous and dissipative hydrodynamic evolution reduces them 
to (8y) = 2.09, 2. 16 and 2.26 at RHIC and 2.82, 2.86 and 2.92 at LHC, respectively. The dissipative 
hydrodynamic results are plotted along with the experimental data from AGS, SPS and RHIC in 
Fig. [| The simple linear extrapolation of the results from the lower energy collisions deviates 
from the RHIC data, which implies that collisions become increasingly transparent as the energy 
becomes higher. The deviation, on the other hand, would be smaller when the hydrodynamic 
effects are taken into account. This indicates that more energy is available for the production of a 
hot medium than that implied from the experimental data. 

Finally, the effects of the cross terms are investigated. As mentioned earlier, the linear scalar 
terms are combined into the effective bulk viscosity; 

n = -Cbn^V M «" - Cna^i + + 
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Figure 2: Mean rapidity losses at E917 (AGS) |g6|, NA49 (SPS) [g7p and BRAHMS (RHIC) [g5[ |28p 
experiments and those of initial distributions and dissipative hydrodynamic results at RHIC and LHC. The 
dotted line denotes a linear extrapolation of the AGS and SPS results and the bars on the RHIC data show 
ambiguities due to the lack of experimental data at forward rapidity [E5[] . 



T T \de 



"BO 



T \dtiBoJ 



(3.2) 



If one assumes £ n<5 „ B = [n B o/(eo + Po)] x CnSe, C = (Cnn + CnSe^)/ T follows from thermody- 
namic relations. It should be noted that the response of a system to expansion, i.e., the bare bulk 
viscosity £nn could be as large as the shear viscosity r\ , but cancellation by the cross term would 
occur except for the crossover region where the sound velocity is small. Phenomenologically, the 
bare bulk viscous term discourages the expansion but the cross term encourages it in an effort to 
decrease the temperature/chemical potential of the system, leading to the cancellation of the effects. 
This could be a reason for the general smallness of the overall bulk viscosity. 

The baryon dissipation current, a vector dissipative current, has the thermo-diffusion term as a 
cross term. Diffusion induced by thermal gradients is known as Soret effect in statistical mechanics. 
This could be potentially as large as the bare baryon dissipation. The numerical result shows, on the 
other hand, the effect is rather small for high-energy heavy ion collisions (Fig. |3|). This is because 
the cross conductivity should vanish at the vanishing limit of the chemical potential because of the 
matter-antimatter symmetry condition V(— jUg) = — V(A*fl)> limiting Soret effect at forward rapidity 
where the net baryon number is non-vanishing. The shear viscosity does not have the cross term 
because it is the only tensor dissipative current. 

The second order cross terms could also be important because the shear stress tensor is gen- 
erally much larger than the bulk pressure, which in turn is larger than the baryon dissipation 
current. It is numerically found that the bulk-shear coupling term in IT and the shear-baryon 
and the bulk-baryon coupling terms in V 11 would potentially be large and should be carefully 
treated for the second order theory to be valid. For example, Fig. |3| (b) shows the effect of the 
baryon-shear term ^y 7r A^ v V p 7r V p in Eq. ( |2T6| ) when the transport coefficient is parametrized as 



Xy K = Cy^riBo/ieo +Pq)]tv with c Vk = 5, and -5. In this case the second order cross term ex- 
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Figure 3: Net baryon rapidity distributions of the color glass-based initial state (solid lines) and of dissipative 
hydrodynamic results with (a) the linear cross factor c v w = 5 (dotted lines), c v w = (dashed lines) and 
c v w — —5 (dash-dotted lines) and (b) the second order cross factor c vlz = 5 (dotted lines), c v % — (dashed 
lines) and c VK — —5 (dash-dotted lines) at RHIC. 



hibits the trend similar to the linear one. Note that the magnitude of the effects are dependent on 
the transport coefficients which are in general not well known. 

4. Discussion and conclusions 

In this study, I developed a relativistic dissipative hydrodynamic model at finite baryon density 
to estimate the net baryon rapidity distribution. The transverse dynamics is integrated out as it is 



motivated by the experiments Q12Q. Numerical estimations are performed with RHIC and LHC 
settings with the equation of state constructed from lattice QCD calculations in Taylor expansion 
method. The initial conditions for the energy and the net baryon number distributions are based on 
the color glass picture. The results suggest that the net baryon number is carried to forward rapidity 
with the flow during the time-evolution, effectively reducing the mean rapidity loss at both RHIC 
and LHC. This could explain the rather high transparency of the collisions at RHIC compared 
with that implied from the simple extrapolation of the AGS and SPS results. They also suggest 
that more kinetic energy would be converted into the inertial energy of the QGP than previously 
deduced from the experimental data. The results at RHIC show that the effects of baryon diffusion 
and viscosity could be visible. The thermo-diffusion effect, or Soret effect, was investigated and 
found to be limited in the forward rapidity region in high-energy heavy ion collisions due to the 
matter-antimatter symmetry condition on the cross conductivity. 

Introduction of transverse dynamics and more realistic transport coefficients would increase 
accuracy of the model. The former would not directly affect the net baryon distribution but it can 
affect the space-time evolution itself and thus the lifetime of the produced hot medium. The latter 
would be important because the dissipative effects could be underestimated without the realistic 
chemical potential dependences of the transport coefficients. It would also be worth-investigating 
to see the correspondence between the mean rapidity losses for the net baryon number analyzed in 
this study and the ones for the energy which has recently measured at LHC [p9|]. 
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